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Abstract: Optical parametric down-conversion (PDC) is a central tool 
in quantum optics experiments. The number of collected down-converted 
modes greatly affects the quality of the produced photon state. We use 
Silicon Photomultiplier (SiPM) number-resolving detectors in order to 
observe the photon-number distribution of a PDC source, and show its 
dependence on the number of collected modes. Additionally, we show 
how the stimulated emission of photons and the partition of photons into 
several modes determine the overall photon number. We present a novel 
analytical model for the optical crosstalk effect in SiPM detectors, and use 
it to analyze the results. 
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1. Introduction 

For more than two decades, the optical process of parametric down-conversion (PDC) HI has 
been the primary photon source in quantum optics experiments IS). This nonlinear process 
converts single photons into pairs of signal and idler photons, which possess quantum correla- 
tions between their various degrees of freedom. The statistical nature of this process is a subject 
of great interest, as the photon-number statistics of the process provide a direct indication of 
its non-classility fSj-jS). Additionally, the photon statistics have practical implications on state 
preparation schemes, where the quality of the produced state depends on the number of photons 
generated simultaneously within the PDC process Q . 

Recent developments in photon-number resolving detectors ll8 UT2ll enable the implementa- 



tion of a wide range of applications in quantum optics, including state preparation schemes 
lUSHTSl, quantum state filtering (161 and quantum gates 117111811 . Number discriminating capa- 
bilities have also been shown to aid in probing eavesdropping attacks on quantum cryptography 
protocols [191. Furthermore, one of the most prominent advantages of number-resolving de- 
tectors is the abiUty to directly measure the photon-number statistics of an optical state. 

The photon-number distribution of a PDC source can range between a thermal distribution, 
when a single down-converted mode is collected, and a Poisson distribution when an infinitely 
large number of spatial and spectral modes (multimode) is collected ll20l . The photon statistics 
of PDC sources has been directly measured in several works ||3]|4l|2T]|22l, all of which report 
a Poisson distribution. However, the majority of quantum optics schemes take great care to 
ascertain that only a single mode is collected. Using an indirect homodyne detection method, 
a thermal distribution from a parametric amplifier was demonstrated 1231 . Transition between 
thermal to Poisson distributions was indirectly observed through the g'^' parameter f24\ and 
through high-order PDC events li25l . In this work, we use a photon-number resolving detector 
in order to directly measure the photon statistics of the signal and idler photons collected from 
a controlled number of modes, allowing us to observe the range of distributions between the 
single-mode and the multimode extremes. 

The photon-number resolving detector used in this work is the Silicon PhotomultipUer 
(SiPM) 1261 . The SiPM detector is composed of multiple avalanche-photo-diodes (APD) op- 
erating in Geiger-mode, combined on a single substrate. When a photon impinges on an APD 
element, an electric discharge can be generated. The photons impinging on the detector are 
distributed across its many elements. The output signals generated by all APDs are combined 
to form a single output pulse whose intensity is proportional to the number of triggered ele- 
ments. In comparison to other technologies, the SiPM detector offers a better photon-number 
resolution (up to 20 photons, compared to 10 or less), easy integration in optical setups and 
operation at room temperature. On the other hand, its photon detection efficiency is 2-6 times 
lower, depending on the working wavelength. 

The photon statistics measured by SiPM detectors may deviate from that of the original state. 
The three main reasons for these deviations are: 1 ) the non-perfect overall photon detection effi- 
ciency (loss), which is determined by the detection and coupling efficiencies, 2) false detections 
caused by thermal excitations (dark counts), and 3) optical crosstalk, in which secondary pho- 
tons, created by carrier relaxation in one APD element, are detected by a neighboring element, 
falsely increasing the photon count l27l . Another deviation is caused by events when two or 
more photons impinge on the same element and are detected as one. However, we have shown 
using a numerical model that this effect becomes significant only when the number of detected 
photons is higher than 20% of the total number of elements 1281 . The descriptions of loss 
and dark counts are quite straightforward. However, several currently available models for the 
crosstalk effect are either computationally difficult l28l or only applicable under limiting con- 
ditions 12914311 . We introduce here an analytical crosstalk model which can be applied for the 
entire range of experimental settings. 

This paper is organized as follows: In Sec.|2]we present a model for the distortion effects in 
the SiPM detector, that includes a new approach to the crosstalk effect. This model is used for 
the interpretation of experimental data and the reconstruction of the original photon-number 
statistics. Sec. [3] describes the experimental setup for the generation of PDC states and their 
detection using an SiPM photon-number resolving detector In Sec.|4]we present measurements 
of photon-number distributions of parametric down-conversion and their dependence on the 
number of collected modes. Conclusions are presented in Sec.|5] 



2. The SiPM detection model 



In order to interpret the photon-number statistics measured using SiPM detectors, the effects 
of loss, dark counts, and optical crosstalk must be properly modeled. If the measured photon- 
number distribution is represented by a vector of probabilities p„„ the original probability dis- 
tribution po, is related to the measured distribution through the relation 

p,„^M-po. (1) 

Each component of the distortion matrix M is the conditional probability M„,„ = M{n\m) of 
measuring n photons given that m original photons arrived at the detector. Using Eq. [T] the 
original distribution can be reconstructed, in principle, by multiplying the measured distribution 
with the inverse of the distortion matrix. This method, however, is limited to relatively large 
probabilities and may produce non-physical results otherwise 132] . Another approach, which 
we use here, is to apply a fitting algorithm which finds the original distribution based on a 
goodness of fit test. 

The matrix M can be written as a product of three matrices, M = M^, • M^/a • M/m v, which 
represent the loss, dark counts (dk) and crosstalk (cf) effects. The individual matrices are con- 
structed as follows: 

Loss 

The probabiUty for loss is described by the relation ll32l 

M/„,,(«|m) = <^ \ n J ' " - (2) 

[O n> m, 

where rj is the overall detection efficiency. 

Dark counts 

The probability of dark counts is assumed to follow a Poisson distribution 1321 of the form 

{0 n < m 

A;r"-exp(-A,,) (3) 
(n— m)! — ' 

where X^k is the average number of dark counts. 
Crosstalk 

Optical crosstalk occurs when a secondary photon generated during an electronic discharge 
in one APD element is detected in a neighboring APD element, thereby creating a secondary 
avalanche. Photons produced in the secondary avalanches can then continue to trigger further 
avalanches in additional elements (crosstalk-generated-crosstalk). The crosstalk probabilities in 
our model are calculated by dividing the crosstalk events into stages and recursively counting 
the number of possible crosstalk events in each stage. The recursion formula for the crosstalk 



matrix element is: 
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where £ is defined as the overall probability for a crosstalk-avalanche to be generated in any 
of the four neighboring elements. Each stage has m triggered elements which can trigger up 
to m new crosstalk events. The probability of generating net new crosstalk events is composed 
of the probability that n^, elements out of the m possible will generate crosstalk, while the 
remaining m — riet elements will not. This probability also includes the combinatorial factor of 
choosing the rict crosstalk-generating elements out of the possible m. The itc crosstalk-triggered 
elements continue to generate additional crosstalk events in a recursive way, until all of the 
n — m crosstalk elements are triggered. 

This novel analytical crosstalk model does not limit the number of crosstalk events and also 
accounts for crosstalk-generated-crosstaUc. The model was developed since other available an- 
alytical models ||29]43TI were inconsistent with our experimental results due to limitations on 
the crosstalk probability values and the number of crosstalk events. Our model does not ac- 
count for the geometrical arrangement of the APD elements in the detector, nor its finite size. 
However, we show below that for values within the experimental range, this analytical model 
agrees with a numerical model that has no such limitations 1 1281 . Additionally, this model is less 
computationally demanding than the numerical model. Thus, it is a useful tool in the interpre- 
tation of experimental data from photon-number resolving detectors, where crosstalk between 
the detection elements plays a crucial role. 

3. The experimental setup 

We have measured the photon-number distributions of the polarization modes of a stimulated 
type-II collinear PDC process. The experimental setup is shown in Fig.[TJa). Pulses of 150 fs 
at a repetition rate of 250 kHz are amplified, frequency doubled to 390 nm and used to pump a 
2 mm long type-II collinear j3 — BaB204 (BBO) nonlinear crystal. The degenerate wavelength 
signal and idler photons at 780 nm, which are created with orthogonal polarizations (horizontal 
(H) and vertical (V)), are split at a polarizing beam splitter (PBS) and detected separately using 
two photon-number resolving detectors. This work focuses on measuring the photon-number 
distribution of a single polarization mode. Thus, the detection of both polarization modes is only 
used to evaluate the heralded detection efficiency of the SiPM detector The results presented 
for a single polarization mode also apply to the photon pair statistics, as a collinear process 
produces the same statistics for pairs and for individual photons. Before the photons are coupled 
to the detectors, they are spectrally and spatially filtered by bandpass filters and by optical fibers, 
respectively. The number of spatial modes which are collected can be changed by using optical 
fibers with different core diameters and numerical apertures ll33l . and the number of collected 
spectral modes can be changed by using filters of different bandwidths. 

The detection configuration is shown in Fig.[Tlb). The optical signal from the photon source 
is coupled to an SiPM detector with 100 elements {Hamamatsu Photonics, S10362-11-100U) 
using optical fibers. The output electrical signal from the detector is amplified using low-noise 
amplifiers, digitized and then analyzed by FPGA electronics, where the number of detected 
photons is extracted in real time. This data is continuously transmitted to a computer, which 
presents the statistics. In order to minimize the effects of dark counts and after-pulsing 1271 . 
the 1 ns long sampling time of the analog electrical signal is synchronized with the arrival time 
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Fig. 1. (a) The experimental setup. A type-II collinear BBO crystal is pumped by 390 nm 
amplified pulses at a repetition rate of 250 kHz. The signal and idler photons are split at a 
Polarizing Beam Splitter (PBS) according to their polarization and detected by two SiPM 
detectors. The number of spatial and spectral modes which are collected is varied using 
interference filters (IF) with different bandwidths and different optical fibers, (b) The de- 
tection configuration. An SiPM detector produces a signal whose intensity is proportional 
to the number of impinging photons. This signal is amplified, digitized using analog-to- 
digital converters (ADC) and analyzed using FPGA electronics in real-time. The data is 
continuously transmitted to a computer, which displays the photon-number distribution. 



of the photons. The dark counts are further reduced by moderately cooling the detectors using 
a thermoelectric cooler to w — 10°C, and the bias voltage level applied to the SiPM is modified 
accordingly. 

4. Results and discussion 

Figure |2] shows an example of a histogram of pulse intensities produced by a coherent state. The 
probability of detecting n photons is proportional to the area of the «''' Gaussian. We obtain high 
photon-number discrimination for as high as 20 photons. The photon-number resolution error, 
which results from the overlapping between neighboring Gaussians, is smaller than 1% below 
12 photons and approaches 12% for 20 photons. The number of resolvable photons is limited 
due to the decrease of the signal to noise ratio as the number of photons is increased. 

We demonstrate the use of our model for reconstructing the original photon statistics with the 
photon-number distribution measurements of a thermal state. The thermal state was generated 
by collecting a single polarization mode of PDC using the setup shown in Fig.[T] The original 
thermal distribution can be written as 

where n is the average number of photons. The measured distribution exhibits a lower value of 
n due to losses and deviates from the power-law dependence of Eq. |5]due to optical crosstalk 
and dark counts. We measured the same thermal state with three different sets of values for the 
amount of loss, dark counts and crosstalk, obtained by operating the SiPM detector with dif- 
ferent bias voltage values ll27l . As a result, we obtained three different measured distributions. 




Intensity of output signal [A.U.] 

Fig. 2. A histogram of the electrical output signal level for a coherent input state. The data 
was accumulated over a period of 30 seconds. Good peak separation is maintained up to 20 
photons. 

although the original photon-number distribution was the same. The measurements are shown 
in Fig. [3] The main distortions in the distributions are the low values of n due to losses and the 
deviation from a straight line in the semi-log plot when the number of measured photons ex- 
ceeds 1 . The latter is the effect of the crosstalk process, which falsely increases the probability 
of measuring high numbers of photons. The crosstalk probabilities of these measurements range 
from the highest probability value (maximum bias voltage), which creates the highest deviation, 
to the minimal probability (minimal bias voltage) which creates the lowest deviation. 

The original distribution is reconstructed by fitting the measured data p,„ to the function 
Pm — Mt(M^//.M/ojj/97/j, where p,/, is the thermal distribution of Eq. |5]and « is the only free 
parameter. The overall detection efficiency Tj, the average number of dark counts A^/^,, and 
the crosstalk probability e, that are used for the construction of the distortion matrices, are 
measured separately. The value of rj is defined as the heralded efficiency between the two down- 
converted polarization modes and is determined by the ratio between two-photon coincidence 
and single-photon counts in the limit n ^ 1 . The dark counts and crosstalk parameters are 
obtained through a separate measurement of the Poissonian dark count statistics which are fit 
to the function p„, = M^M^itl, where 1 is a vector of zeros with 1 at the first position, and X^k 
and e are the free fit parameters. 

All three measurements in Fig. [3] were reconstructed to the same original state, within the 
margins of error, despite the different measured distributions. This result indicates the consis- 
tency of our evaluation of the dark counts, loss and crosstalk parameters. We compare these 
results to those obtained using the numerical crosstalk model introduced in Ref. Il28l . The 
numerical model is based on simulations of the detection process, and takes into account the 
geometrical configuration of the detector, the finite number of detection elements, and pos- 
sible attempts to trigger neighboring cells which have already been triggered. The model is 
characterized through a parameter e„„, defined as the probability of generating crosstalk in one 
particular neighbor among the four nearest neighbors, rather than the overall probabiUty of 
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Fig. 3. Measurements of a thermal state conducted on a polarization mode of a type- 
II collinear PDC source, spatially and spectrally filtered using a single-mode fiber and 
a 3 nm bandpass filter. The experimental data was fitted to the obtained statistics using 
the model presented in this work (solid lines). The reconstruction of the original pho- 
ton statistics from the measured data is also shown. The data was obtained at a tem- 
perature of — 10°C. The bias voltage values and the fit parameters are: Vbias = 68.3V, 
T] = 6 • 10"^ ± I • 10"^ Arf^ = 2.4 ■ 10^3 ± 6 ■ 10-'', 7j = 8.9 ± 0.5, e = 0.280 ± 0.005 (green 
circles), V^ia^ = 67.7 V.T] = 4.9- 10"^ ±3- 10-^, = i.Q. lO^^ ± 3 • lO"'*, 7? = 8.8±0.5, 
e = 0.140 ±0.006 (red upward triangles), and Vtias = 61 AN , rj = 4.2 - 10"^ ±5 • lO""*, 
l^k = 2.5- 10"'' ± I • 10"'', M = 7.8 ±2.7, e = 0.040 ±0.006 (blue downward triangles). 
The corresponding crosstalk values obtained with the numerical model of Ref. 1281 are 
£„„ = 0.078±0.001, e„„ = 0.038±0.001, and e„„ = 0.0I0±0.001, respectively. Fits using 
the numerical model are presented by black dashed lines. 

generating crosstalk among the nearest neighbors, which we define here as e. The values of e,„, 
obtained from the fit correspond to the values of e through the relation £ = 1 — (1 — 1281 . 
The agreement between the two models shows that although the analytical model does not ac- 
count for the geometry of the detector, it provides a good description of the crosstalk process 
for all crosstalk values in the experimental range. 

We now consider the photon-number distributions of down-converted photons as a func- 
tion of the number of collected modes. If « photons are distributed among s modes, such that 
ni +n2 + ■ ■ .ris = n, then the overall probability of measuring « photons is given by 1201 : 



where the product nj^j/:'(«,) is the joint probability of measuring n,- photons in each mode i, 
and the combinatorial factor amounts to the number of possible arrangements of n indistin- 
guishable photons into s modes. If the photons are distributed evenly among the modes, the 
average number of photons in each mode is uniform and equals «7 = n/i. The photon-number 
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Fig. 4. (a) Representative phioton-number distribution measurements and their fits. The re- 
sults are presented for distributions which range between thermal and Poissonian statis- 
tics. Experimental errors are smaller than their respective symbol size. The average pho- 
ton number for each data set is reduced as less modes are collected, (b-e) A close-up of 
each of the measured distributions along with distributions calculated for the same fit pa- 
rameters while changing the mode-number value between a single, 1.4 and 5.5 modes, 
and a Poisson distribution. The measurements were taken with the following parameters: 
ri = 0.013, Xdk = 1 .95 • 10"^, e = 0.26, 3 nm bandwidth filter, single mode fiber for 780 nm, 
integration time of 780 minutes (green circles), rj = 0.019, X^f, = 2.1 • 10^^, e = 0.23, 
3nm bandwidth filter, single mode fiber for 1500 nm, integration time of 100 minutes 
(blue squares), rj = 0.013, Xjk = 1-69 • 10"^ e = 0.27, 10 nm bandwidth filter, single 
mode fiber for 1500 nm, integration time of 100 minutes (orange triangles), and rj = 0.020, 
^dk = 3 • 10^^, £ = 0.22, 3nm bandwidth filter, multimode fiber for 780 nm, integration 
time of 115 minutes (red diamonds). 



distribution within each of the modes is thermal, /:>(«,) = Pth{ni), and by substituting Eq.|5]into 
Eq. |6] we obtain the following probability distribution 

("'t"7')(TTl)W' ™ 

known as the negative binomial distribution. This distribution converges to a Poisson distribu- 
tion when the number of modes approaches infinity. 

The dependence of the photon-number distribution on the number of modes is demonstrated 
by measuring the photon statistics collected from a controlled number of modes. Figure|4]shows 
four representative measurements of the photon-number distributions which range between a 
single mode (thermal) distribution, and a multimode (Poisson) distribution. The solid lines are 
fits to the function p,n = MctM^i^MiussPs, where ps is defined in Eq.|7]with n and s as the free 
parameters. We detected up to 14 photons and obtained good fits to the data. Notice that the 
difference between the distributions becomes more marked as the number of photons increases 
(see Figs. |4tb)-(e)). The difference is almost undetectable for low photon numbers. In fact, 
measurements of at least 10 photons are required for our experimental parameter range, in 
order to properly discriminate between distributions with different mode numbers. 

The number of modes obtained for all combinations of fibers and filters is shown in ta- 
ble [T] The non-integer mode numbers should be interpreted as a weighted number of occupied 



Table 1. The number of modes as obtained from the photon-number statistics. The table 
shows the number of modes for interference tilters with different bandwidths and for optical 
fibers with different mode field diameters (MFD) and numerical apertures (NA). Measure- 
ments in the first row were taken with a single-mode fiber for 780 nm, the second row with 
a single-mode fiber for 1550 nm and the third row with a standard graded-index multimode 
fiber. Distributions are considered Poissonian for s> 10. 
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modes and result from a non-uniform distribution of photons between the modes. We can dis- 
criminate between distributions containing up to about 10 modes, as distributions containing 
over 10 modes cannot be distinguished from Poisson distributions. Table [T] shows the pro- 
portional increase of the number of spectral modes with the increase in the bandwidth of the 
bandpass filters. The width of the phase-matched spectrum of the down-converted photons is 
estimated to be « lOnm. Due to spectral distinguishability, this width includes about 3 temporal 
modes Il24l|25ll34l . The number of spatial modes measured with our system is also consistent 
with the choice of fibers. When the PDC photons, which passed through a 3 nm bandpass fil- 
ter, were coupled into a single mode fiber for 780 nm, a thermal distribution was produced. 
Whereas, when the photons were coupled into a multimode fiber, a Poisson distribution was 
produced, even though the same filter was used. The relatively small number of modes ob- 
tained with the single-mode fiber for 1550 nm can be attributed to non-uniform mode coupling, 
which does not equally stimulate all of the modes in the fiber. 

The stimulated nature of PDC can also be used to obtain the number of collected modes 
ll25l . The average number of down-converted photons in a single mode is determined by the 
stimulation parameter T,, which has a linear dependence on the pump field, the crystal length 
and the nonlinear coefficient of the crystal ll35l . Assuming an equal distribution of the photons 
between s modes, the mean number of photons n is given by 

n ~ s -Tii 

= i • sinh^(T,) 

= s ■ sinh- l^aVfj , (8) 

where Wi is the average number of photons in each of the individual modes, a is the coupling 
constant between the pump field and the nonlinear crystal, and / is the pump intensity. 

In order to show the dependence of « on the number of modes, we measured photon-number 
distributions as a function of the pump intensity for the different numbers of collected modes. 
The values of « were extracted by fitting the different distributions to Eq. [T] The fit values 
of n as a function of the total pump power are shown in Fig. |5] All measurements show a 
nonlinear dependence, indicative of a stimulated process. The number of modes s obtained from 
these fits corresponds well to the number of modes obtained through fits to the photon-number 
distributions (see Fig.|4|. Furthermore, all fits resulted in similar values of the a parameter, in 
agreement with the model of Eq.[8] 
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Fig. 5. The average number of photons collected from multiple modes as a function 
of the pump intensity. The solid lines are fits to Eq. [8] with the following parameters: 
(a) s = 1.1 it 0.7, a = 0.08 ±0.02 (green circles) (b) s = 1.4 ±0.5, a = 0.10 ±0.02 
(blue squares) (c) s = 63 ± 0.9, a = 0.08 ± 0.01 (orange triangles), (d) .s = 23 ± 8, 
a = 0.06 ± 0.01 (red diamonds). 



5. Conclusions 

In conclusion, we have built a detection setup, which incorporates a Silicon Photomultiplier 
(SiPM) as a photon-number resolving detector Our setup enables good photon-number dis- 
crimination for as many as 20 photons. The output signal of the SiPM is analyzed in real time 
using FPGA electronics. The measured photon-number distribution probabilities are interpreted 
by modeling the effects of loss, dark counts and optical crosstalk, which distort the original 
photon-number probabilities. We present a useful analytical crosstalk model that is applicable 
for all crosstalk probability values in the experimental range and includes crosstalk-generated- 
crosstalk. 

We have performed measurements of photon-number distributions of a single polarization 
mode of a collinear type-II PDC process. The high photon-number resolution of our system en- 
ables the differentiation between distributions collected from a different number of spectral and 
spatial modes. We present distributions which range from single-mode (thermal) to multimode 
(Poissonian) statistics. The results show how the photon-number distribution is determined by 
the number of collected modes. Furthermore, we show the dependence of the average number 
of photons on the number of collected modes, and the stimulated nature of parametric down- 
conversion. 



